THE GROUND STATE ENERGY OF THE THREE DIMENSIONAL 
GINZBURG-LANDAU MODEL IN THE MIXED PHASE 

AYMAN KACHMAR A B 

Abstract. We consider the Ginzburg-Landau functional defined over a bounded and smooth 
three dimensional domain. Supposing that the strength of the applied magnetic field varies 
between the first and second critical fields, in such a way that Hex ^ H <C Hc 2 , we estimate 
the ground state energy to leading order as the Ginzburg-Landau parameter tends to infinity. 



1. Introduction and main result 

We consider a bounded and open set C I 3 with smooth boundary. We suppose that f2 
models a superconducting sample subject to an applied external magnetic field. The energy of 
the sample is given by the Ginzburg-Landau functional, 

S 3D (iP, A) = ty, A) = J (\(V- MA)^f + y(l - H 2 ) 2 ) dx 

+ k 2 H 2 I |curlA-/3| 2 (ia;. (1.1) 

Here k and H are two positive parameters, the wave function (order parameter) tp £ // 1 (f2;C), 
the induced magnetic potential A 6 i^ iv F (]R 3 ), where i^ ivF (M 3 ) is the space introduced in 
(|1.2p below. Finally, f3 is the external magnetic field that we choose constant, (3 = (0,0, 1). 

Let if 1 (M 3 ) be the homogeneous sobolev space, i.e. the closure of C^°(1R 3 ) under the norm 
u ^ ll u ll_ffi(R3) := ll^ tt lli 2 (R 3 )- Let further F(x) = (— X2/2, xi/2, 0). Clearly divF = 0. 

We define the space, 

^div,F( M3 ) = { A : divA = 0, and A - F e i^QR 3 )} . (1.2) 
Critical points (tp, A) £ H 2 (Q; C) x Hi, F (M 3 ) of £ 3D satisfy the Ginzburg-Landau equations, 
' -(V - inHA) 2 ifj = k 2 (1 - |V>| 2 )V> in Q 

< curl 2 A = - — Im^V - inH A)^)l n in M 3 (1.3) 

, v ■ (V - inHA.)ip = on <9s2 , 

where In is the characteristic function of the domain 0, and is the pointing interior unit normal 
vector of d£l. 

For a solution (ijj, A) of (|1.3p . the function ^ describes the superconducting properties of the 
material and //curl A gives the induced magnetic field. The number k is a material parameter, 
and the number H is the intensity of a constant magnetic field externally applied to the sample. 

In the mathematics literature, Type II superconductors usually correspond to the limit k — > 00, 
see |H [18] . In this regime one distinguishes three critical values Hc\ , Hq 2 and Hc 3 for the applied 
field. Those critical fields are roughly described as follows. If H < Hc x , the material is in the 
superconducting phase. Mathematically, this corresponds to > for any minimizer (^,A) 
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of (jl.ip . If Hc 1 < H < Hq 2 , the magnetic field penetrates the sample in quantized vortices 
(corresponding to zeros of tp). If Hc 2 < H < Hq 3 , superconductivity is confined to the surface of 
the sample (corresponding to \ip\ very small in the bulk). Finally, if H > Hc 3 , superconductivity 
is lost, which is reflected by ip = everywhere in £1. In this paper, we will focus on the regime 
when the applied magnetic field varies between Hc\ an d Hq 2 . In the scaling we choose in this 
paper, this regime corresponds to Iiik/h <C H <C k as k — > oo. Here, if a(n) and 6(k) are two 
positive functions, the notation a(n) <C b(n) means that a{n)/b{K) — > as k — > oo. 

In the case of two dimensional domains, which correspond to infinite cylindrical superconduct- 
ing samples, there exists a quite satisfactory analysis of the critical fields i?Cu Hq 2 and Hc 3 - 
As we can not give an exhaustive list of references, we invite the reader to see the monographs 
[H [18], where a detailed review of the material is present. Still in the two dimensional setting, 
the most accurate available characterization of the critical field Hc 2 is given in |13^ I14| . 

The situation is less understood in three dimensions, especially the regime of magnetic fields 
close to the first critical field Hc t - For a superconductor occupying a ball domain, a candidate 
for the expression of the critical field Hq 1 is given in [2|. Related results are obtained for 
superconducting shells in [7]. For general domains, the analysis of the critical field Hq 3 started 
in |16j . then a sharp characterization of Hc 3 is given in |1U| . In the papers [HE], it is proved that 
superconductivity is confined to the surface of the domain, provided that magnetic field is close 
to and below Hq 3 - A fine characterization of the critical field Hc 2 together with leading order 
estimates of the ground state energy in large magnetic fields are recently obtained in |12| I15j . 
This paper is complementary to those in |12^ [T5] , 

The ground state energy of the functional in (]l.ip is defined as follows, 

C (k,H) =inf{£: 3D ^,A) : A) G H\n;C) x H^R 3 )} . (1.4) 

The main result of this paper is Theorem II . II b elow . It is a generalization of an analogous result 
proved for the two-dimensional functional in |20j. 



Theorem 1.1. Suppose that the magnetic field H is a function of n and satisfies 

hiK 

C ii C as k — y oo . 

K 

Then, the ground state energy in (|1.4|) satisfies, 

Co(k, H) = \Q\ kH In + o(kH In \J~^^j > as K ~ * 00 • (1-5) 

As immediate consequences of Theorem 11.11 we obtain that, if (if), A) is a minimizer of (|1.1|) . 
then the induced magnetic field curl A is close to the applied magnetic field /3, and that the mag- 
nitude of the order parameter \^\ is close to 1 almost everywhere in Q. The physical meaning of 
this is that the applied magnetic field penetrates the sample almost everywhere and concentrates 
along 'vortex lines'. On these vortex lines the order parameter ift is expected to have zeros (this 
is not rigorously proved in this paper), but away of them, the sample remains in the supercon- 
ducting phase (|^| is close to 1). Therefore, the regime considered in Theorem 11.11 corresponds 
to what is actually named in the physics literature as the mixed phase. 

In the course of the proof of Theorem 11.11 we obtain the following conclusions as immediate 
corollaries. 
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Corollary 1.2. Under the assumptions made in Theorem \l.ll if (if), A) £ H (Q;C) x i/j iv (R 3 ) 
is a minimizer of the energy in (jl.ip . then, as k — > oo 7 




-Here, dx zs i/ie Lebesgue measure in Q, the measure e K ^n{$>, A) and i/ie current A) are 

respecizve/y, 



In two dimensions, /i Kj _f/ is a measure and it is proved that it gives the density of vortices, 
hence it is called the vorticity measure, see |18) . 

The proof of Theorem 11.11 is obtained as follows. First we start by the analysis of an approx- 
imate problem in a 'large' cube. The cube geometry allows us to link this problem to another 
two dimensional problem in a square. The later is analyzed using tools from |18] , 

Using a ground state of the approximate problem, we construct a test configuration whose 
energy provides an upper bound of the ground state energy Cq(k,H). As a consequence of this 
upper bound, we obtain that, for a minimizer (ip, A) of (|1.1|) . the induced magnetic field curl A 
is close to the applied field f3 in L 2 -norm. Using this and the regularity of the curl-div system 
in R 3 , we get an estimate of A — F in C 0,1 / 2 -norm. 

The a priori estimates obtained for minimizers allow us to determine a lower bound of the 
energy that matches with the obtained upper bound. Actually, we use the 'semi-classical' local- 
ization techniques developed in [8] to reduce the problem to that of the approximate problem in 
a cube. Then the analysis of the later problem is used to obtain the matching lower bound. 

An interesting aspect of the analysis is that we do not use constructions involving vortices, i.e. 
we do not localize the set where {i£l! : ^(x)! < 1/2} (as this is certainly difficult in three 
dimensions). This is a significant difference between the strategy of our proof and the one given 
in |20] for the two-dimensional functional. However, the construction of 'vortex-balls' for the two 
dimensional functional 'implicitly' appears in the analysis of the three dimensional approximate 
problem, as we refer to results of |18[ I20| . In the context of the Ginzburg-Landau model, the 
implementation of 'semi-classical' techniques to address situations where vortices exist seems 
rather new. 

The analysis presented in Section [2] combined with a recently proved estimate in |12] enables us 
to prove a theorem of independent interest (Theorem 12.41 below) . which concerns the asymptotic 
behavior of a limiting constant appearing in |19j . thereby answering a question raised by the 
authors of the aforementioned paper. 

The paper is organized as follows. Section [2] is devoted to the analysis of the approximate 
problem. In Section [3j an upper bound of the ground state energy is obtained. In Section [4j 
interesting estimates are obtained for minimizers of (jl.ip . Section [5] is devoted to the proof of 
the lower bound. 
Remark on notation: 

• The letter C denotes a positive constant that is independent of the parameters k and H, 
and whose value may change from line to line. 
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• If a(re) and are two functions with b(n) 7^ 0, we write a(n) ~ if o(k)/6(k) — )• 1 
as k — y 00. 



2. The approximate problem 

2.1. Two dimensional energy. Let K = (—1/2,1/2) x (—1/2,1/2) be a square of unit side 
length, h cx and e be two positive parameters. Consider the functional defined for all u G 

tf^-c), 

£ 2D (u) = jf_ jj(V-i/i cx A )u| 2 + ^(l-M 2 ) 2 ) dx. (2.1) 
Here Ao is the vector potential, 

A (xi,x 2 ) = ^(-x 2 ,xi) , (xi,x 2 )GR 2 , (2.2) 

whose curl is equal to 1. 

Notice that the functional E 2T) is a simplified version of the full Ginzburg-Landau functional 
considered in [20j , as the magnetic potential in (|2.ip is given and not an unknown of the problem. 

We introduce the ground state energy, 

m (h ex ,e) = mi{E 2D {u) : u G H l {K;C)} . (2.3) 

Since E 2T) is bounded from below, there exists a ground state (minimizer) associated to 
m o(h>ex, e)- If u is such a ground state, then it results from a standard application of the maximum 
principle that, 

|u| < 1 in K . (2.4) 

Consider the regime of magnetic fields h ex as in Theorem 12.11 below. We can obtain a lower 
bound of mo(h cx ,e) (or rather of E 2Y) (u), with u a ground state) exactly as in |18| Section 8.2], 
by using a scaling argument that reduces the situation to magnetic fields of lower order (precisely 
of order | lne|). In this way, we get the following theorem. 

Theorem 2.1. Assume that h ex is a function of e such that 

I In el /i cx <C —z , as e — > . 

Then the ground state energy mo(/i ex ,e) satisfies, 

m (h cx ,e) > h ex \n—==(l + o(l)) , 
ev n ex 

as e — > 0. 

Minimization of the functional E 2I) over 'magnetic periodic' functions appears naturally as 
well. Let us introduce the following space, 

E h ^ = {ueHl c (R 2 ;C) : u( Xl + 1, x 2 ) = e ih ^ 2 u{x u x 2 ) , 

u( Xl ,x 2 + 1) = e- ih ™ Xl/2 u{ Xl , x 2 )} , (2.5) 

together with the ground state energy, 

m p (/i cx ,e) = mf{E 2D (u) : u G E he J . (2.6) 
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Theorem 2.2. Assume that /i ex is a function of e such that 

I In el <C /i cx <C -o , as e — )■ . 
TTien //ie ground state energy m p (h ex ,e) satisfies, 

m p (h ex ,£) = /lex hi — t==(1 + o(1)) , 

eV/lex 

as e — >■ 0. 

Proof. Since the restriction of a function in E^ to A is a function in ff 1 (A'), we get that 
m p {h cx ,e) > mo(h ex , e), where mo(h ex ,e) is the ground state energy in (|2.3p . Theorem 12. II then 
gives us a lower bound of m p (h ex ,e). 

We prove the upper bound by computing the energy of a test function u constructed in [6]. 
Let N be the largest positive integer satisfying N < ^//i ex /2vr < N + 1. Divide the square iiT 
into iV 2 disjoint squares (-Kj)o<j<jv 2 -i each of side length equal to 1/N and center aj. 

Let h be the unique solution of the problem, 

—Ah + /i ex = 27r5 ao in Kq 
dh 

^ = on 5A 

ay 

/ /idx = 0. 

Here v is the unit outward normal vector of Kq. By uniqueness of h as solution of the afore- 
mentioned problem, h is symmetric with respect to the axes of the square Kq and hence satisfies 
periodic conditions on the boundary of Kq. Moreover, the function v(x) = h(x) — ln\x — ao\ is 
smooth in Kq, since — Av + /i ex = 0. Consequently, through a scaling argument, it is easy to 
check that, as e — > 0, 



/ \Vh\ 2 dx < 2^1n-^- + 0(1) 

JK \B(ao,e) 



< 2^ In — l — + 0{l). 

£ y /l ex 



We extend h by periodicity in the square K. Let <f> be a function (defined modulo 27r) satisfying 
in K \{aj : < j < N 2 - 1}, 

V<p = -V ± h + /lexA • 
Here = (— d X2 ,d Xl ) and Ao is the magnetic potential in (|2.2p . 

If a; € -f^O) let p(x) = min(l, |x — ao|/e). We extend the function p by periodicity in the square 
K. We put u{x) = p{x)e^ for all x£ K. Then ii can be extended as a function in the space 
Eh cx hi (|2.5p . see [5| Lemma 5.11] for details. 

The energy of it is easily computed, since u is 'magnetic periodic'. Actually, 



E 2D (u) = N 2 x J k (p 2 |V/i| 2 + |Vp| 2 + ^(1 - p) 2 ) dx 
<N 2 x ( 2vr In — ^= + 0(1) ] . 



Since, AT = y/h^j2TT^l + o(l)) as e — > 0, and m p (h cx ,e) < £/ 2D (u), we deduce that, 

m p (h ex ,e) < /i cx In —^=(1 + o(l)) , 



as e -+ 0. □ 
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Since mo(h cx ,e) < m p (h ex ,e), we get as a corollary of Theorems 12,11 and 12.21 

Corollary 2.3. Let mo(h ex ,e) be the ground state energy introduced in (|2,3p above. Suppose 
that h ex is a function of e and | lne| <C h ex <C 1/e 2 as e — > 0. Then, 

m (h cx ,e) = h cx ln — L=(l + o(l)) , 

as e — >• 0. 

Theorem 12,21 serves in answering a question of independent interest arising in |19] . Consider 
two constants b G (0, 1) and R > 0. Let K R = (-R/2, R/2) x (-12/2, R/2). If u G H l (K R ), we 
define the energy, 



'K R 

together with the ground state energy 



Fk r {u) = ^ (b\V - iA )u\ 2 + \u\ 2 )^j dx . 



e p (b, R) = M{F Kr (u) : u G E R ] . 

Here Ao is the magnetic potential introduced in (|2.2p and E R is the space introduced in (|2.5|) . 
(with h ex = R). It is proved that, for all b G (0, 1), there exists a constant f(b) such that, 

/<»)-5«m*M (2.7) 

The limiting constant f(b) appeared in |19|. IT] . then it is recently studied with different tools in 
|12j . This limiting constant describes the ground state energy of both two and three dimensional 
superconductors subject to high magnetic fields (see [12J). 

The behavior of the function f(b) as b — > 1_ is analyzed in details in |12j . However, the 
behavior as b — > 0+ remains open. Only a non-optimal estimate on f(b) is given as b — > 0_|_ in 
|19| . Here, using Theorem 12.21 and an estimate in |12j . we describe the leading order asymptotic 
behavior of f(b) as b — > 0_)_. 

Theorem 2.4. Let f(b) be as defined in (|2.7p . Then, as b — > 0+, fib) satisfies, 

f{b)= b - In -L (1 + (1)). 

Remark 2.5. In |19J . it is proved that 

b-^<f(b)< b -ln±= (l + o(l)), 



as b — > 0_i_ . 



Proof of Theorem \2.4\ It is proved in |12} Theorem 2.1 &; Proposition 2.8] that there exist uni- 
versal constants C and i?o such that, 



V 6 G (0,1), V R > Rq 



2/(6)- e ^ R) 



(2-8) 



R? 

Let h cx = R 2 and e = y/b/R. A scaling argument shows that, 

e p (b,R) = bm p (h cx ,e). (2.9) 

We select R = 1/6 so that as 6 — > 0+ we have e — > and | lne| <C h cx <S e~ 2 . Theorem 12.21 then 
tells us that, 

m p (h cx ,e) = h cx ln—== (l +o(l)) . 
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We insert this estimate into (|2.9p then we substitute the values h ex = R 2 and £\/h ex = \fb. 
Finally, inserting the resulting estimate into (|2.8p finishes the proof of the proposition. □ 

2.2. Three dimensional energy. If T> is an open set of M 3 and u G H^D; C), we define the 
energy 

G v (u) = j (b\{V -i¥)u\ 2 + ^{l-\u\ 2 ) 2 ^j dx. (2.10) 

Here F is the magnetic potential, 

F(xi, x 2 , £3) = (— X2/2, xi/2, 0) , { Xl ,x 2 ,x 3 ) £K 3 , (2.11) 

whose curl is equal to 1. 

Let 6 and R be two positive parameters. Consider a cube Qr of side length R defined as 
follows, 

Q R = (-R/2, R/2) x (-R/2, R/2) x (-12/2, R/2) . (2.12) 

We introduce the ground state energy, 

M (b,R) =mi{G QR (u) : u G H\Q R ; C)} . (2.13) 

In the next theorem, we give an asymptotic lower bound of the ground state energy Mq (6,12) as 
b — > and R — > 00 simultaneously, in such a way that IniRb' 1 / 2 ) < 12 2 . 

Theorem 2.6. Suppose that the positive parameters b = 6(e) and R = 12(e) are functions of a 
parameter e such that, 

6(e) ^0, 12(e) ^00, and — ^ In -> , 

12(e) 2 yjb{e) 

as e — > 0. 

Then, the ground state energy Mo(b,R) satisfies, 

M (b,R) 1 , 

— ^— = 61n -^(l + o(l)), 

as e — )• 0. 

Proof. Let /i cx = 12 2 and e = \fb~jR. By the assumption on 6 and 12, it is easy to see that e — > 
and I lne| <C /i cx < 1/e 2 . 

Consequently, Corollary 12.11 tells us that the ground state energy mo(h ex ,e) in (|2.3p satisfies, 

m (h DX ,e) = /i cx ln — 7==(l + o(l)) . 

We will prove that, 

M (b,R) = bRm (h cx ,£), 

which will immediately give us the asymptotic estimate in Theorem 12.61 

Let u G H l {Q R -C), K = (-1/2,1/2) x (-1/2,1/2) and Q 1 = K x (-1/2,1/2). Define the 
rescaled function u G H (Qx; C) as follows, 

V x G Qi , 5(x) = u{Rx) . 

It is easy to check that, 

G Qr (u) = 612 (|(V - i/i cx )F)5| 2 + ^ (1 - |5| 2 ) 2 ) dx ± ) dx 3 



-1/2 

-1/2 



> 



612 y ^ (jf (KV^ - i/i ex )F)5| 2 + ^ (1 - |5| 2 ) 2 ) dx ± j dx 3 . 
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Here, if x = (rci, x%, x%) G M 3 , we write xj_ = (x\,X2) and V x± = (d xl ,d X2 ). Then, recalling the 
definition of mo(h ex ,e), we get, 

Gq r (u) >bR m (h ex ,e)dx 3 = bRm (h cx ,e) . 

J -1/2 

Taking the infimum over all functions u G H^{C^yi',^S), we get that i\/fo(6, R) ^ bRTTtQ(h cx ,s). 
Let Uh^e be a ground state of E 2T) , i.e. £' 2D (n/i exi£ ) = mo(h ex ,e). Define the function, 

u:Q R 3x^ u hcxj£ (x ± /R) . 

Then, G Qj! (ti) = bRE 2D (u hc ^ e ), thereby showing that M (b,R) < bRm (h cx ,e). □ 

3. Upper bound of the energy 
The aim of this section is to give an upper bound on the ground state energy Cq(k,H) in 

dm). 

Theorem 3.1. Assume that the magnetic field H satisfies Iiik/k <C H <C k as k ^ oo. T/ien 
i/je ground state energy Cq(k,H) in (jl.4|) satisfies, 

C {k, H) < \U\kH In (l + o(l)), (3.1) 

as k — ^ oo. 

Furthermore, there exists a constant kq such that, if n > kq and (ip, A) is a minimizer of the 
functional in (jl.ip . then 



||curl(A-F)|| i2(R3) < ^LJhx^. (3.2) 

Proof. Notice that if (ip, A) is a minimizer of HUD, then £ 3D (V>, A) = C {k,H). Consequently, 
the estimate in (|3.2p follows immediately from the upper bound in (|3.ip . 



/ K H\ 1/4 1 

Let b = H/k and £ = ( ) , Then, as k — > oo, we have, 

Vin^y 

k<l, £ < 1 , ^\/^H : > 1 . 

Let /i ex = l/^ 2 and e = y/bi. Then, as k — > oo, we have £<1 and | lne| <C h cx <C l/e J . 

Recall the ground state energy m p (h cx ,e) and the space introduced in (|2.6p and (|2.5 

respectively. Let u G -Eft cx be a ground state corresponding to m p (/i ex ,e), i.e. 



J M(V - ^exA )n| 2 + - |u| 2 ) 2 ) drr 



m p (h ex ,e) 



For all x = (xj_,X3) G M 3 , we introduce the function, 

v(x) = u(iV kH x±) 
Let (Qj) be a lattice of ffi 3 generated by the cube, 

. 1 1 n / 1 i W l 



2£Vk;# 2£\/kHJ V 2£Vk# 2W kH J \ 2N kH 2WkH 
It is easy to check that, 



\{V - inHF)v\ 2 + ^-{1 - \v\ 2 ) 2 ) dx = —^=m p {h cx ,e). 
* J WkH 
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Here F is the magnetic potential in (j2TTj) . Let J = {Qj : Qj n d£l / 0} and N = Card J. 
Then, as k — > oo, we have, 

AT = |0| x (ly/^H) 3 (1 + o(l)) . 

Recall the functional £ 3D in (jl.ip . We compute the energy of the test configuration (t>, F). Since 
curlF = f3 and the function v is magnetic periodic with respect to the lattice Qj, we get, 

£ 3D (v,F) = Nx J (j(V - i K HF)v\ 2 + y(l - \v\ 2 ) 2 ^j dx 

= N x - 7 L=?7i p (/i cx ,e) . 
We use Theorem 12.21 the definitions of h ex and e, and the asymptotic behavior of TV" to get, 

N x _L=m p (/i e x,£) = Kj ffln JlL (l + o(l)) , 
1\/kH V 

as k — ^ oo. This proves the upper bound of Theorem 13.11 □ 



4. A PRIORI ESTIMATES OF MINIMIZERS 

The aim of this section is to give a priori estimates on the solutions of the Ginzburg-Landau 
equations (jl.3p . Those estimates play an essential role in controlling the error resulting from 
various approximations. 

The starting point is the following L°°-bound resulting from the maximum principle. Actually, 
if (V>, A) G iJ^fijC) x if] ivF (IR 3 ) is a solution of ([Eg]) , then 



U»(0)<1- (4-1) 
Next we prove an estimate on the induced magnetic potential. 

Proposition 4.1. Suppose that the magnetic field H is a function of k such that In k <C kH <C k 2 

as k — y oo. There exist positive constants kq and C such that, if k> kq and (ip, A) S i? 1 (f2; C) x 
if^ ivF (]R 3 ) is a minimizer of the energy in (jl.ip . then, 



A - F\\ H 2 {n) < Jin 



Vtf V V H 



||A-F|| c0 , 1/2( n)<^yiny^. 
Here F is the magnetic potential introduced in (|2,lip . 

Proof. The estimate in C 0,1 / 2 -norm is a consequence of the Sobolev embedding of H 2 (Q) in 
Notice that it follows from Theorem 13.11 that. 



curl(A-F)|| L2(R 3) < -&Wln^, ||(V - miTA)^|| L2(n) < 2|n|v^ff Win ^ . (4.2) 



Let a = A — F. We will prove that ||a||#2(Q) < -^===^/ln \j — ■ Since diva = 0, we get by 



regularity of the curl-div system (see e.g. [U Theorem D.3.1]), 

l|a||L6(R3) < C|| curl a|| L 2 (R3) . (4.3) 
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The second equation in (|1.3p reads as follows, 

- Aa = — Imhp (V - inH AU)l n . 
kH 

Select a positive constant M such that the open ball K = B(0,M) contains Q. By elliptic 
estimates (see e.g. [3J Theorem E.4.2]), 

ll a llH 2 (fi) < C(II«IIl 2 (x) + ll^o|| L 2( X )) . 

Using the embedding of L 2 (K) into L 6 (K), the estimate in (|4.3p and the bound < 1, we get 
that, 

I|a||j9*(n) < C (\\ curla|| £a(R s) + ^IK V ~ « k - h ')^IIl 2 (^)) • 
Inserting the estimates in (|4.2|) into this upper bound finishes the proof of the proposition. □ 

5. Lower bound of the energy 

In this section, we suppose that D is an open set with smooth boundary such that D C 0. 
We will give a lower bound of the energy, 

£ (1>,A;D) = ^|(V-miiA)^| 2 + y(l-^| 2 ) 2 ) dx , (5.1) 

where (ip,A) is a minimizer of the functional in (|l.ip . The precise statement is the subject of 
the next theorem. 

Theorem 5.1. Suppose that the magnetic field H is a function of k such that Inn <C kH <C k 2 
as k — y oo. If (ifi, A) € i? 1 (0;C) x i^ iv F (IR 3 ) is a minimizer of the function in (|l.ip . then, 

£^,A;D) > \D\kH In J^ + oLhIu J^) , 



as k — y oo. Here £o(ip, A; D) is introduced in (|5.ip . 

Proof. Let ^ € (0, 1) be a parameter (depending on k) that will be chosen later in such a way 
that (^/kH)~ 1 «f«l as k — y oo. Consider a lattice (Qj)j of ]R 3 generated by the cube, 

Q e = (-1/2,1/2) x (-1/2,1/2) x (-1/2,1/2) . 

Let i7 = { j : Qj C I?} and N = Cardj7\ Then, as k — > oo, the natural number N satisfies, 

N = W + °{l>)' (5 - 2) 

Moreover, we have the lower bound, 

fb^AjDj^^fo^AjQj). (5.3) 

For each j £ JT, we will bound from below the term £o(ip, A; QA Let Xj be the center of the 
cube Qj. Using the estimate of ||A — F|| c o,i/2/q) given in Proposition 14. 1| we may write for all 
x E Qj, 

\A(x) - F(x) - (A( Xj ) - F(xj))\ < C\i 1 ' 2 , 
where C is a constant that is independent of j, x and n, and the parameter A is defined by, 
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We define <fij(x) = (A(xj) — F(xj)) ■ x, uj(x) = e % ^ x 'ip(x) and a.j(x) = A(x) — Vipj(x). Then 
we may write, 

VieQj, |a,-(x) -F(x)| < CX£ 1/2 , (5.5) 

and 

So(il),A;Q j )=S (u j ,aj;Q j ). (5.6) 

We may write, for all 5 € (0, 1), 



|(V - iKHaj) Uj \ 2 > (1 - <S)|(V - we£TF)u,-| 2 - 2cT 1 (K#) 2 |a i 



We insert this estimate into the expression of £o(uj,a.j; Qj) then we use the estimate in (|5.5 
and that \uj\ = to get, 



£ («i,aj-; Qj) > (1 - *)£b(uj,F; Qj) - C5~ 1 {kH) 2 \ 2 1 I \^\ 
1 kH and b = H/k. For all i£K 3 such that \x\ < R, we define, 



dx . 



(5.7) 



Let R = 

vj (x) = u I X 
Then a simple change of variable shows that, 

£o(uj,F; Qj) -- 



x 



1 



„ G Q*( v j)> 
by kH 

where Gq r is the functional in (|2.10|) and Qr is the cube in (|2.12p . 
We select £ in the following way, 



(5.8) 



In k J 



(5.9) 

<C 1 as k — > oo. Conse- 



1 R 

With this choice, we have (\/ kH)" 1 < I < 1, 1 < and — ^ In — = 

V ^ 

quently, Theorem 12.61 tells us that the ground state Mo(b,R) in (|2.13p satisfies 

M (6, fl) = 6fi 3 In —= (1 + o(l)) . 
Vb 

Since «j 6 H 1 (Qji), we get GQ R (vj) > Mo(b,R). Substituting this into (|5.8|) and using the 
aforementioned asymptotic expansion of Mo(b,R), we get, 



£b(«i,F;G i ) = -^In^=(l + o(l)). 
By inserting (|5.10p into (|5.7p and using (|5.6|) . we get for all j € 

A; Qj-) > (1 - <*)-_ In -= (l + o(l)) - C(T^(«F) 2 A 2 / |# 



(5.10) 



Taking the sum over j ^ J and using 



we get, 



£(V>, A; D) > (1 - S) N x -^L In (l + o(l)) - Crt(ftif) 2 A 2 / 



dx . 



(5.11) 



where N = CardJ". To finish the proof, we use the bound < 1, the definition of A in (|5.4|) . and 
we choose 5 = I 1 / 2 . This gives that the remainder term in (|5.1ip is equal to o(nH In yj 'k/H ) . For 
the leading order term in (|5.1ip , we use the asymptotic expansion of N in (|5.2p , that R = l\f kH , 
and we observe that it is equal to 



effing (1 + (1)). 
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□ 

Proof of Theorem \l.l[ Combining the upper bound in Theorem 13.11 and the lower bound in 
Theorem 15.11 with D = Q, we get the estimate of the ground state energy in Theorem 11.11 □ 

Proof of Corollary The convergence of curl A — ft in L 2 (]R 3 ;R 3 ) is proved in Theorem 13.11 
Since div(A - F) = and A - F G ij 1 (R 3 ; M 3 ), we get that, 



< 

we infer from Theorem I3JJ 



||Vcurl(a-F)|| L 2 (R 3) = || curl(A - F)|| L2( 

Consequently, it results from the convergence of curl A in I? (M 3 ) that curl A -> /3 in H 1 ( 

We prove the convergence of ^^(t/jjA). Let B(x) = curlA(x). Since divA = 0, it results 
by taking the curl on both sides of the second equation in (ll.3p . 

- AB + B = ii Ki H(i>, A) in . 

Since B -»• /3 in i/ 1 (]R 3 ; R 3 ), we get that -AB + B -> (3 dx in £>'(]R 3 ;R 3 ). 

It remains to prove the convergence of the measure e^jf^, A). It suffices to prove that 
e Ki ij(V>, A) — > dx in the sense of measures. If D is any open set in Q with smooth boundary, 
then we have by Theorem 15.11 . 

£ {^,A;D) > \D\kH In (l + o(l)) , £o(ip, A;Q\D) > \U \ D\kH In (l + o(l)) . 

Here £q(4>, A;D) is introduced in (|5.ip . Recall the functional £ 3D in (jl.ip . Since 
£ (iP, A; D) + Soty, A;Q\D) = £ 3D (%b, A) 

KH\n\hiJ^(l + o(l)), 

1 

£ (iP, A; D) = \D\kH hxJ^(l + o(l)) . 
This is sufficient to conclude the convergence of e Ki i^('i/ ; , A) to in the sense of measures. □ 
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